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Abstract 

Nonlinear integral equations are proposed for the description of the full 
finite size spectrum of the 2-dimensional 0(3) nonlinear cr-model in a periodic 
box. Numerical results for the energy eigenvalues are compared to the rotator 
spectrum and perturbation theory for small volumes and with the recently 
proposed generalized Liischer formulas at large volumes. 



1 Introduction 



The study of finite size (FS) effects in quantum field theories has recently been re- 
ceived a lot of attention. While understanding the structure of FS effects has always 
been an important part of the theory of quantum fields (in particular in the numeri- 
cal simulation of lattice field theories), this renewed interest is due to the important 
role FS effects are playing in the verification of the AdS/CFT correspondence pp. 
Partially motivated by similarity to the AdS / CFT problem, nonlinear integral equa- 
tions (NLIE) have been proposed [2] to describe the spectrum of excited states in the 
2-dimensional 0(4) nonlinear a- model confined to a finite, periodic box. Motivated 
also by a problem in the AdS/CFT correspondence, a generalization of Luscher's 
formulas [3j, giving the leading large volume dependence for all excited states were 
proposed [4]. This result was succsessfully applied for calculating the 4- loop [4] and 
5-loop [5] anomalous dimension of an important operator in the AdS/CFT corre- 
spondence. 

In this paper we propose a set of nonlinear integral equations and corresponding 
quantization conditions for the description of the full finite size spectrum of the 0(3) 
NLS model. Based on the proposal that the 0(3) NLS model can be represented as 
a limit of appropriately perturbed parafermion conformal field theories (CFT), 
the Thermodynamic Bethe Ansatz (TBA) equations for the ground state of the 
0(3) NLS model were proposed firscl in [7]. The equations could be formulated 
in a rather elegant form in terms of the incidence matrix Ij k of an infinite Dynkin 
diagram of type V depicted in Figure 1: 

oo 

\ogy k {e) = -£5 k0 coake + ^2l kj (K*hyYj)(6), Y k {6) = l + y k (6), k = 0,1,2,... 

3=0 

(1.D 

where i is the volumcl, K(6) denotes the TBA kernel K(8) = 27r cosh(fl) anc ^ * denotes 

oo 

convolution: (/ * g)(x) = f dy f(x — y)g(y). The energy of the ground state can 



1 Later a generalization to a more general class of G/H coset models were given in [6]. 
2 We measure all energies and lengths in units defined by the infinite volume mass gap of the 
model. 
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Figure 1: Dynkin-diagram associated with the Y-system of the 0(3) cr-model. 
be expressed in terms of the solution of (11. ip by a simple integral expression: 

oo 

E = — [ — coshtf \nY (6). (1.2) 

y 27r 

— oo 

It was also shown that the solutions of the ground state TBA equations satisfy cer- 
tain functional equations, the so-called Y-system equations, which take the following 
form for the 0(3) NLS model: 

oo 

y k (6 + iq)y k {e-iq) = X[Y 3 {6) I ^ ) g=| k = 0, 1, .... (1.3) 

3=0 

Detailed studies of analogous problems for integrable lattice models [8] indicated 
that the Y-system is universal in the sense that different solutions of the same Y- 
system describe all excited states of the model. The difference between the solutions 
is encoded into their analytic structure. Thus the excited states can also be described 
by TBA equations, which are similar to the ground state TBA equations but are 
supplemented by additional source terms and quantization conditions. 

For relativistic quantum field theories (QFT) this phenomenon was first observed 
in [9], where for certain perturbed CFTs it was demonstrated that excited state 
TBA equations can be obtained by the analytic continuation of some parameters 
of the ground state equations. Later, based on an integrable lattice regularization, 
universality of the Y-system for all excited states was demonstrated and excited 
state TBA equations were derived in the sine-Gordon model [10]. In this model it 
was also shown that the source terms for excitations can be determined by studying 
the large volume asymptotics of the Y-functions. 

The above examples suggested that it is true in general also for relativistic QFTs 
that the same Y-system describes all excited states of the model and the TBA 
equations for the ground state and excited states differ only by source terms and 
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additional quantization conditions, which can be determined from the knowledge of 
the infrared (IR) asymptotics of the Y-functions. Accepting the above hypothesis 
1-particle TBA equations were proposed for the 0(3) and 0(4) NLS models by an 
appropriate modification of the analytic properties of the vacuum Y-functions [llj. 

The disadvantage of the TBA approach to sigma models is that the number of 
unknown functions is infinite. It is desirable to find an equivalent nonlinear integral 
equation (NLIE) description of the problem, with only a few unknown functions. 
Although such a set of NLIEs is equivalent to the Y-system and TBA, it is a very 
useful technical tool that makes the analysis of finite size effects more efficient. In 
the cases of interest a Y-system is equivalent to a T-system and the latter can be 
solved by an auxiliary linear problem called TQ-relations. As it was shown in a 
series of papers the TQ-relations are the basic tool for the derivation of an NLIE 
[12j . In [13} [HJ US] integrable lattice regularizations were used to derive NLIEs for 
certain QFTs, while in the case of the 0(4) NLS model [2] started directly from the 
TQ-relations of the continuum theory. 

In a series of integrable lattice models p2] and relativistic QFTs [13j [HI [15] it 
has also been demonstrated that similarly to the TBA case the NLIE for the ground 
state and excited states differ only in source terms and quantization conditions, 
which can be determined from the IR analysis of the unknown functions. Using this 
hypothesis and starting from the ground state equations [THl [TT] 1-particle NLIEs 
were proposed for the 0(3) and 0(4) NLS models [T8] . 

The organization of the paper is as follows. In the next section we propose a set of 
nonlinear integral equations and quantization conditions for the full spectrum of the 
periodic 0(3) NLS model. In Section 3 we study the large volume limit and give a 
detailed analysis of the spectrum of 2-particle states in this limit using the solution 
of the asymptotic Bethe Ansatz. In Section 4 we discuss our numerical results, 
both for small volumes, where comparison to the rotator spectrum and perturbation 
theory, and for large volumes, where comparison with the Bajnok-Janik formulas, 
are presented. Technical details are explained in Appendices A and B. 
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2 The nonlinear integral equations 



Our starting point in this paper is the hypothesis that the same Y-system describes 
all excited states of the 0(3) NLS model. We generalize the NLIE technique to 
describe the complete finite size spectrum of the model. The generalization is based 
on the equivalent formulations of the Y-system (jl.3p through a T-system and TQ- 



relations. With the help of the TQ-relations we can define appropriate auxiliary 
functions in terms of which the infinite set of TBA equations can be reduced to a 
finite set of equations, containing only a few unknown functions. 

In the large I limit eigenstates of the model are multi-particle states satisfying the 
Bethe-Yang equations, which provide quantization conditions on the set of particle 
rapidities, X = {9^}, k = 1,2, ... ,N. This is constructed as follows. We take an 
M-magnon solution of the asymptotic Bethe equations: 



M, (2.1) 



ipo(uj + iiT) _Q (uj + i7t) ._. 
<p (uj -m) Qo(uj - i 7r) ' 

where 

N M 

mo) = - Ok), Q o (0) = - %)> ( 2 - 2 ) 

fe=l j=l 

and the UjS are the magnon rapidities and their number M < N determines the 
isospin of the iV-particle state through the formula J = N — M. 
Next we construct the T-system elements for k — 0, 1, . . . 

k 

T k {9) = Qo(6 + i(k + l)q)Q o (0 - i(k + l)q) ]T £(0 + i(k - 2j)q), (2.3) 

j=0 



where 



^ {) Q (e + ig)Q (e - i q y [ ' ) 



In particular, 

T o (0) = MO + n) MO - iq), T x {6) = <p (6) f x {9), (2.5) 

where 

f ^°) = 7T7ff\ (MO + wr) Qo(0 - m) + MO - m) Qo(0 + m)} , (2.6) 
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which is a polynomial of degree N as a consequence of the Bethe Ansatz equations 
( 12.11) . It can be shown also that all T k are polynomials of degree 2N. 
Using the definition 

<l>(e) = (p (e)(p (e + 2iq) (2.7) 

it is possible to verify that the functions <fi, To, Ti, . . . form a semi-infinite T-system 

T k (9 + i q) T k {9-iq)- T fc _i(0) T fc+1 (0) = <p(9 + ikq) <j>(0 -ikq), (2.8) 

where is the complex conjug of 0. The system is semi-infinite because it can 
be consistently truncated putting T_i{6) = 0. 
Next one constructs the Y-system elements 

y k (9) = Tk-i(Q)T k+ i{6) ^ ^ 

n + ikq) (j){x — ikq)' 



Y k {6) = T fc( g + ^) T fc( g ~^) ; ( 2 1Q ) 

n + ikq) <p(x — ikq)' 



which satisfy the Y-system equations 

Vk(0 + iq) y k (0 - iq) = Y k+1 (9) Y k ^(9) (2.11) 

for k = 2, 3, . . . and 

yi (e + iq) yi (9-iq)=Y 2 (e). (2.12) 
The main result in this construction is that 

yi (9) = A(6 + iq\6), (2.13) 

where A(6,6_) denotes the eigenvalue of the transfer matrix made out of the two- 
particle S-matrix Sij(9) of the 0(3) NLS model: 

T(9,9) = Tr o (£oi(0 - #i) S 02 (9 - 9 2 )...S 0N (9 - 9 N )). (2.14) 

More precisely, here we deal with the spin-1 transfer matrix, where both in the 
auxiliary space and in the "quantum" spaces we take the spin-1 representation of 
SU(2). 



^Throughout this paper / denotes the complex conjugate of the function /. 



5 



It is a common experience in the theory of Bethe Ansatz equations that the 
polynomials T k have only real roots inside the physical strip \lm8\ < q (see, for 
example, [H]). We denote the set of real roots of T k (which may be empty) by t k 
for k — 1, 2, From (I2.5P we see that 

ti = Xufi, (2.15) 

where t± is the set of real roots of T\ and from (12. 9p we see that, denoting the set of 
roots of y k inside the physical strip by r/ k , 

% = ViUt w , fc = 2,3, ... Vi=h- (2.16) 

Furthermore, we see that the quantization conditions 

Y k (z - iq) = 0, zet k , k = 2, 3, . . . Y x {z - iq) = 0, z G t~i (2.17) 

are satisfied by (12.101) . 

We need one more equation to quantize the rapidities of the physical particles. 
They are determined from yo(9) by the equation: 

Y o {0 h -iq) = Q, k = l,2,...N (2.18) 

For large i this equation is equivalent to the Bethe- Yang equations 

e iis ' mMk A(9 k \9) = -1, (2.19) 

if yo(9) = e~ £coshe 2/1(6*). This relation is consistent with the Y-system equations fll.3p 
independently of the value of I. In this paper we will assume that this relation holds 
exactly. 

For finite i the Y-system equations (12.111) for k = 3,4, .. . and (I2.12p are un- 
changed and only the k = 2 equation is modified to 

V2(e + iq) y 2 (6 - iq) = Y 3 (9) Y x {9) Y (6). (2.20) 

Since the only change is the inclusion of the multiplyer Yq, which is very close to 
unity for large £, we will assume that the finite I solution is a smooth deformation 
of the Bethe Ansatz construction, i.e. both the functions y k (6) and the sets t k are 
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smooth deformations of the ones determined by the above explicit construction at 
large t. (12.1 5p still holds, but with X containing the modified particle rapidities 9 k . 

Now it would be possible to rewrite the Y-system equations for excited states as a 
set of TBA integral equations similar to (II. ip . by adding source terms corresponding 
to the yk roots. This would be supplemented by quantization conditions of the form 
( I2.17P and (I2.18p . In this paper we leave out this step and proceed directly to the 
NLIE description. 

Assuming that the solution of the excited state Y-system (II. 3p is found, we first 
build a T-system, which is equivalent to it. Since the semi-infinte Y-system (12. lip 
for k — 3, 4, . . . is of the form of the standard A\ Y-system, we can find suitable T- 

system functions 0, T k satisfying (12.81) so that (12.91) and (12.101) hold for k = 2, 3, 

The k = 2 relation (12.201) corresponds to the modified relations 

Y 1 (9) Y (9) - 1 = yi (9) [Y (9) + e~ ic ° she ] = — g^§^ , (2.21) 

-■>'" + %q) <p{x — iq) 



^W = S^. (2-22) 
+ iq) (p{x — iq) 



In Appendix A it is shown explicitly that it is always possible to find T-system 
functions T k for k = 0, 1, . . . such that ft2~2T\f . (12T221 are satisfied, ([23]) and fT2TT0]) 
hold for k = 2,3... and the roots of T k coincide with the elements of for k = 
1,2,... Furthermore the T-system equations (I2.8P are satisfied for k = 1, 2, . . . with 
(f) given by (12. 7p (with the modified 9^ rapidities). The T-system functions obtained 
this way are smooth deformations of the corresponding ones constructed in the large 
I limit. 

Next we construct the TQ system. This is based on the fact that the T-system 
equations (I2.8P are integrable and have a Lax representation [2] through the auxiliary 
problem (TQ-relations) 

T k+1 (9)Q(9+ikq)-T k (9-iq)Q(9+i(k + 2)q)=<j)(9+ikq)Q(9-i(k + 2)q), 
T k ^ 1 (9)Q(9-i(k+2)q) - T k (9-iq)Q(9-tkq) = -0(9-tkq)Q(9+ikq). (2.23) 

It can be shown that in the Bethe Ansatz limit Q(9) = Q(9) = Qq(9) satisfies 
the TQ system (I2.23p . In Appendix A we show how to find a Q{9) solution which 
is a smooth deformation of Qo(9). 
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Both equations (12.81) and (I2.23P are invariant under the gauge transformation 

T k {6) -> g(9 + ikq)g(9-ikq)T k (9), 
m g(e + iq)g{9-iq)<f>{0), 

Q(Q) _ g(e-iq)Q{6). (2.24) 

The Y-functions y k are gauge invariant. 

We will use for the NLIE the following gauge invariant auxiliary functions: 

Q(9 + l (k + 2)q)T k (9-tq) 

bk{6) = Q(e-i(k + 2) q )<K9 + ik q y "- ' (2 ' 25) 
BfcW = g(e-i(fc + 2) g )0((9 + i]fc g ) fc - ' (2 ' 26) 

where 

B fc (0) = l + 6 fc (0). (2.27) 
Furthermore the auxiliary functions (I2.25f2.26p satisfy the functional equations 

bk(9)h(9) = Y k (0)[l + 8 kl y o (9)], k>l, (2.28) 
B k (9 + iq)B k (9-iq) = Y k+1 (9), k > 1. (2.29) 

(I2.27p . (12.281) and (12.291) . together with the analytic properties that can be read off 
the reperesentation (12. 25112. 26p are the key to find a set of NLIEs that effectively 
truncate the TBA equations at the kth node. 

The complete set of unknown functions of our NLIE we discuss below are the 
functions: b%(9), which will be denoted by b{9 — ij) for short, y\{9) and yo(9). For 
the quantization conditions the function bo(9) will be used as well. 

The derivation of the NLIE is given in Appendix B, for an important subset of 
multiparticle states. Here we present the equations and quantization conditions for 
the most general excited state of the model precisely. The NLIEs take the form: 

log b{9) = iTt5 b + iD + "i(9)+ig+~<(9)+ig+ 1 (9) + (G*ln£)(0) - (G+ 27 * In B)(9) 
+ (K~^ *]nY 1 ){9) + (K~^ +r ' *lny o )(0), S b 6 {0, 1} 

log yi (0) = iirSy + ig y {9) + (A + ^ * \nB)(9) + (A^+ 7 * \nB)(9), 8 y e {0, 1} 
log y (9) = -i cosh 9 + \ogyi(9), 

B(9) = l + b(9), Y 1 (9) = l + y l (9), Y {9) = 1 + y {9). (2.30) 
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where 7 is a contour shifting parameter restricted into the interval < 7 < |-, 
In denotes the "fundamental" logarithm function having its branch cut on the neg- 
ative real axis and we introduced the notation for any function /: 

The kernel function G of (I2.30p reads as 

G < 9 > = / 1 ^Sr (2 - 31) 

—00 

while K{9) is the kernel of the TBA equations. The source terms of (12.301) read as: 

j=i j=i j=i 

M c M w 

~ YtXQ-Cil-YtXiiiP-Ws), (2.32) 

3=1 3=1 

91(8) = ( 2 - 33 ) 

3=1 



( 7r 

g y {6) = \im g y [9 + i- - ir] 



N 2 ^ M s 

UQ) = YlxKie-h^ + J^ixKie-v^ + XKie-v^-J^ixKie-s^ + XKie-s,)) 

3=1 3=1 3=1 

M c M w 

~ Y^Xk{6-c 3 )-Y,Xkii{9-w 3 ), (2.34) 
3=1 3=1 

where x{9) an d Xk{0) are proportional to the odd primitives of the kernel functions 



x (6) = 2tt J dx G(x), Xk(0) = 2vr J dx K(x), (2.35) 


and the second determination of any function fu(9) is defined by: 

fniO) = f{6) + f(6 - 2vrsign(Im^)). (2.36) 
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The function Xk{6) is given by the formula: 

sinh [%- 4- -) 

XK(0)=iln , )* ; , Xk(0) = -xk(-0) WeC. (2.37) 
sinh (z| - |j 

The branch cuts are chosen to run parallel to the real axis so that Xk{9) is an odd 
real analytic function on the entire complex plane and continuous along the real 
axis. In this case Xk(9) is not periodic anymore with respect to 2ni. It is periodic 
only modulo 2tt, i.e. the following identity holds: 

XK (O + 2m)= X K(0)-2n. 

It follows that the distance between the consecutive cuts is 27U and the jump of 
Xk{9) is equal to — 2n at each branch cut crossed from below. The choice of branch 
cuts is depicted in Figure 2. 



5'm/2 


Im 


3'm/2 


in/2 


Re 




-in/2 


-3'm/2 



Figure 2: Locations of the branch cuts of xk(&)- 

Although all physical quantities should be gauge- invariant, i.e. invariant under 
(I2.24p . the description of the source objects is nevertheless more transparent in a 
particular gauge fixed by the condition (12. 7p . In this gauge it is possible to show 
(see Appendix A) that Ti(8j) = for all physical rapidity 8j, j = 1,2, ... ,N and 
we can define T\{9) by 

T 1 {e)=f 1 (e)(p {6). (2.38) 
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In what follows we will work in this fixed gauge. 

The objects appearing in the source terms of (I2.30p are as follows. 

1. Type I holesft : {hf^}, j = 1, . . . ,Ni, which are zeroes of T%(x) satisfying the 
condition: < \lmhj \ < |. 

2. il/o/ej^: {hj}, j = 1,...,N 2 corresponding to zeroes of T 2 (x), with < 
|Im hj\ < 7 

3. "Close objects": {cj} j = 1,...,M C , which are "source objects" satisfying 
the condition: 7 < |ImCj| < 7T. 

4. " Wide objects": {wj}, j = 1, . . . , M w , with 7r < |ImWj|. 

Close and wide objects appear in complex conjugate pairs and they are related 
to the zeroes of the Q function in the same manner as in the spin-1 XXX chains 

pa 123 . 

There are also two types of special objects in the source terms of the equations. 
They are defined by the relations: 

Imbg6-^) = 27r/ i ., \b^{ Sj )\>l, (Imlog6^)'(^) < 0, j = l,...,N s , 

(2.39) 

lmlogb-i(vj) = 2TrI v ., \b^(v j )\>l, (Imlogr^'K) > 0, j = l,...,N$ 

(2.40) 

and they are called ordinary and virtual special objects respectively [20]. The func- 
tion D(Q) of f)2.30p is expressed by the asymptotic rapidities of the iV particle state: 

N 

D(6) = J2xK(0-0 k )- 

fc=i 

The values of the constants of the NLIE (12.301) depend on the number of different 

The term "type I hole" comes from the fact that in the thermodynamic limit of the equations 
the real zeroes of Ti(#) correspond to holes in the distribution of real Bethe roots. 

5 The term " hole" comes from the fact that in the thermodynamic limit of the equations the 
real zeroes of T2(9) correspond to holes in the distribution of 2-strings. 
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source objects: 

5 b = + J (mod 2), and S v = J + M w (mod 2), (2.41) 

where J denotes the isospin of the state. In addition to (I2.30p we need two other 
equations for the determination of type I holes and wide objects. For the determi- 
nation of type I holes, a counting function can be defined: 

Z x {e) = g y {6) + -{K-~<*lnB){6)--{K + ~<*lnB){6)-5 y , < |Im0| < -, (2.42) 

where Zi(6) = —i logb Q (6). The function necessary for the determination of wide 
objects is 

log&o(0 + ig) = igui{9) + (GJ/ * In B)(6) - (G+ 7 *\nB)(9), it < Im^.43) 

Taking into account the fact that Kjj(6) and Xkii{@) are zero, this equation looks 
as if it were the second determination of the first equation in (12.301) . 

The source objects appearing in the NLIE are not arbitrary parameters, they 
have to satisfy certain quantization conditions dictated by the analytic properties 
of the unknown functions of the NLIE. The quantization conditions for "magnonic" 
degrees of freedom read as: 

• For holes: 

\ logft-Tfo) = 2tt4„ j = l,...,iV 2 . (2.44) 

• For close source objects (only for the upper part of the close pair): 

I log&^(cj) = 2nI cU j = l,...,M c /2. (2.45) 

• For wide objects (upper part of the wide pair): 

-i \ogb {w] + iq) = 2nI wU j = 1, M w /2. (2.46) 

(We have determined only the upper parts of the complex pairs, the lower parts are 
simply given by complex conjugation.) 

For ordinary and virtual special objects the defining relations (12.391) and (I2.40|) 
themselves play the role of quantization conditions. 
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• Finally for type I holes: 

Z l (hf ) ) = 2irl h m, j = l,...,N 1 . (2.47) 

In finite volume the rapidities of physical particles are also quantized. Their 
quantized values can be determined from yo(6): 

-logy {6 k + £) = 27rl k , k = l,...,N. (2.48) 
% 2 

All the above quantum numbers I a . 's are half integers. A state is then identified 
by a choice of the quantum numbers (i^., I c ., ...). We note that the NLIE itself 
can impose constraints on the allowed values of the magnonic quantum numbers. 
Moreover there are relations among the numbers of different species of source objects. 
They are called counting equations and in this model there are two of them: 

N 2 + 2N$-2N S = M C + 2J + 2M W , (2.49) 

N 1 -2N S R = J + M 1 - M R , (2.50) 

where Mi denotes the number of close and wide objects lying farther than | from 
the real axis, Mr stands for the number of reajf) Bethe roots of the asymptotic Bethe 
Ansatz equations of the model, and N R means the number of the real special objects. 
Real special objects can be either real Bethe roots or real type I holes. They are 
called specials because at the positions of real specials the counting function of real 
Bethe roots is no more monotonically increasing: i.e. ^-Zi(x) < 0. We note that 
the counting equation (12.501) is valid when all the type I holes are real, while (I2.49P 
is valid also for complex holes as far as their imaginary parts lie within the interval 
[-7,7] 

In the TBA analysis of the sine-Gordon theory [10J it turned out that the en- 
ergy expression for an iV-particle state consists of two terms: the sum of the kinetic 
energies of N non-interacting particles plus an integral expression similar to (11. 2p 
containing the Y-function corresponding to the massive node of the Y-system. As- 
suming this formula to be valid also for the 0(3) NLS model the energy can be 

6 Real roots are those zeroes of 1 + 60 {&) along the real axis which are not type I holes. They 
correspond to real solutions of the asymptotic Bethe equations. 
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expressed in terms of the solutions of the NLIE: 



N r do 

E = J2 cosh0 k- / — cosh 8 \nY (6). (2.51) 

— oo 

3 NLIE in the IR limit and solution of the asymp- 
totic Bethe equations for 2-particle states 

In the I —>■ oo limit the NLIE for N particles reduces to the NLIE of an inhomoge- 
neous N-site spin-1 577(2) vertex model [ISJED], where the roles of inhomogeneities 
are played by the physical rapidities. The NLIE description for the inhomogeneous 
spin-1 SU(2) vertex model is equivalent to the Bethe Ansatz solution, thus in the 
infrared (IR) limit the NLIE (12.301) is equivalent to the asymptotic Bethe equa- 
tions (12. ip . The relation between the magnon rapidities in (12.11) and the source 
objects of the NLIE is as follows. The magnons with |Imwj| > 5 + 7 can be ob- 
tained by the formula Uj = Uj + z | sign(Im [7j), where Uj is the infrared limit of a 
close or wide object. Magnons with |Imttj| < § + 7 are summed up by the integral 
terms of the NLIE this is why no source terms can be associated to them. 

In this limit the quantization condition (I2.48P is equivalent to the asymptotic 
Bethe- Yang quantization for the physical rapidities: 

e u«*e> MOk-in) Q (e k + in) x fc = _ 

ipo{8 k + 171) Qo{0k - ITT) 

The term multiplying e ^ sinh6 ' fe on the left hand side of (13.11) is nothing but the 
eigenvalue of the "color" transfer matrix the 0(3) NLS model taken at the point Q k . 

M d k-irc)Qo(9 k + i7r) 

— 1 a \ ■ \ 7T7a ~: = H^k^)- (3.2) 

Thus the rapidity quantization conditions given by the NLIE ( I2.301l2~48]) in the IR 
limit are identical to the Bethe- Yang equations of the 0(3) NLS model. 



2-particle states 

In this subsection we analyse the structure of the asymptotic solution of the NLIE 
equations based on the explicit T-system solution (12.31) . We consider zero momen- 
tum 2-particle states, with states corresponding to J = 0, 1, 2 in the isospin space. 
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The analytic properties of the problem are determined by the magnonic Bethe roots 
and the zeroes of the transfer matrices T\{9) and T2(9) inside the physical strid_. 
^From (12.31) we see that the above data depend on J and the two physical rapidi- 
ties 8i t 2 and for zero-momentum (symmetric) states with fixed J the state can be 
characterized by the magnitude of 9\. 

The result of the calculation is as follows. 

J = 2 state: 

For all real values of the rapidities 6\ and 9 2 : 

• Qo(0) = 1, 

• T\{9) has two real zeroes, 

• T 2 (9) has four real zeroes. 
J — 1 state: 

For any symmetric state 9\ = —9 2 . Qo(9) = 9, and the distribution of zeroes of the 
transfer matrices depend on the magnitude of 9\. 

< 9 1 < ^tt ~ 2.72070 

• T\{9) has no zeroes, 

• T 2 {9) has two real zeroes. 
^7T < 01 < 7T ~ 3.14159 

• T\{9) has two complex zeroes, 

• T 2 (9) has two real zeroes. 



7T < 9 l < J | + \Jf it ~ 5.78682 



• Ti(9) has two real zeroes, 

• T 2 {9) has two real zeroes. 

7 The strip \Jxa6\ < 5 in the complex plane. 



15 



• Ti(9) has two real zeroes, 

• T 2 {9) has two real and two complex zeroes. 
^pn < 9 1 

• Ti(9) has two real zeroes, 

• T 2 (9) has four real zeroes and there is a special object in the corresponding 
NLIE as well. 

J — state: 

For this state the Qo function in the IR limit reads as: 

Q (9) = (9-u + i()(9 -u - u = (3.3) 

where ( = \\J ^j- + . Depending on the magnitude of the rapidities for this 

state there are seven regions corresponding to the various possibilities for the dis- 
tribution of the relevant zeroes of the transfer matrices T\ and T 2 . In our numerical 
investigations for J = we considered only the simplest symmetric states with 
< 9i < y/ ^-v 7 ^! -^ ~ 1.55809. In this case neither Ti nor T 2 has zeroes in the 
physical strip. 

4 Numerical results 

We have considered the NLIE description of the cases N = (vacuum) and N = 1 
(mass gap) previously [TT], [18]. Some numerical results for these cases are given in 
the Eo, E\ columns of Table [H 

We now discuss the N = 2 (2-particle) cases extensively. This is the next non- 
trivial case, with possible isospin values J = 0,1,2. Concretely, we consider only 
zero momentum states with physical rapidities 9\p = ±H with the smallest possi- 
ble momentum values. In addition, we also consider some special N = 3, J = 1 
3-particle states with physical rapidities #13 = ±H, 9 2 = 0, which leads to an 
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NLIE very similar to the one corresponding to the Rllb region for 2-particle states 
discussed below. 

In all cases we started from the large £ BY solution and found the numerical 
solution of the NLIE equation by iteration. We used the following parameters: 
h = 0.04 (step size for the Simpson formula in rapidity space) 
h' = 0.004 (step size for the Fourier integrals) 
A = 240 (cutoff in rapidity space) 

e = 10- 12 ((absolute) accuracy for solving the quantization conditions) 

With these parameters, it was necessary to perform a few hundred iterations 
(corresponding to less than 1 hour CPU on a laptop). In this way a (relative) 
accuracy of about 2 • 10~ 8 is achieved for both the energy values and the position of 
the roots, except very close to the boundary points between regions where qualitative 
changes occur (discussed below). 

For N = 2, J = 2 there is just one region: starting from large £, we can 
gradually come down to the UV region £ ~ 10~ 6 without any qualitative change in 
the algorithm. The qualitative description is correctly given by the BY-equations: 
there are two pairs of real roots for T 2 (arranged symmetrically around the origin) 
and one pair of real roots for T\ (also symmetrical). Some numerical results are 
given in Table [1] (£22 column) . 

The case N = 2, J = 1 is more complex. From the large £ BY equations we see 
that in this case only T2 has a single pair of real roots. Let us denote this region 
Rlla. From the BY equations we get the estimate £ > 0.74, below which T\ has (a 
complex conjugate pair of) imaginary roots inside the physical strip. Indeed, we can 
monitor the imaginary T\ roots as they move towards the physical strip by solving 
the Rlla NLIE equations for decreasing £ values and we estimate that the boundary 
of Rlla is at £ = 0.65, below which (region Rile) we have to solve the NLIE with 
one complex conjugate pair of imaginary roots for T\ and one pair of real roots for 
T 2 . The limits of this region are 

0.37 <£< 0.65 (4.1) 

(as opposed to the 0.49 < £ < 0.74 estimate from the BY solution). 

At £ = 0.37 the imaginary T\ roots meet at the origin and as £ decreases further 
they move away symmetrically from the origin along the real axis. We denote this 
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£ 


£o 


£i 


£22 


£21 


0.000001 


-0.96553941 


-0.64792555 


-0.01287029 


11.59612929 


0.000003 


-0.96027050 


-0.62282893 


0.05183255 


11.60073719 


0.00001 


-0.95362228 


-0.59128817 


0.13308151 


11.60648275 


0.00003 


-0.94655984 


-0.55793920 


0.21890220 


11.61250075 


0.0001 


-0.93741883 


-0.51502162 


0.32920544 


11.62015568 


0.0003 


-0.92740876 


-0.46835316 


0.44895233 


11.62835934 


0.001 


-0.91395474 


-0.40619591 


0.60808718 


11.63907639 


0.003 


-0.89850538 


-0.33567869 


0.78805971 


11.65090156 


0.01 


-0.87635700 


-0.23640710 


1.04016160 


11.66669458 


0.03 


-0.84838500 


-0.11495754 


1.34580545 


11.68338729 


0.1 


-0.80069857 


0.07772108 


1.82023294 


11.69620904 


1.0 


-0.48624957 


1.08420867 


3.96980668 


11.60661538 



Table 1: Finite size spectrum of the 0(3) model. The Eq — I Eq column corresponds to the 
vacuum state, the £\ — iE\ column to the (standing) one-particle state, and the two last columns 
correspond to the zero-momentum N = 2, J = 2 and N = 2, J = 1 states, respectively. 

region by Rllb. According to the BY solution, Rllb ends at £ = 0.036 but the BY 
equations are no longer relevant for such small £ values and in fact we find that Rllb 
extends all the way to the UV limit £ — > 0. The three regions of the N = 2, J = 1 
problem are summarized in Table [2J Some numerical results are given in Table [T] 
(£21 column). 



region 


range 


T\ roots 


T 2 roots 


Rlla 


£ > 0.65 




1 real pair 


Rile 


0.37 < £< 0.65 


1 imaginary pair 


1 real pair 


Rllb 


£ < 0.37 


1 real pair 


1 real pair 



Table 2: The regions in the N — 2, J = 1 case. 



For N = 2, J = we have not considered all regions. The large £ RlOa region is 
characterized by no roots for either T\ or T 2 . It corresponds to £ > 2, below which 
an imaginary pair of T 2 roots enters the physical strip. The iV = 2, J = problem 
corresponds to a number of different regions, but this is not discussed here in detail. 
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Rotator spectrum 

At small physical volumes, I — > 0, the only relevant degrees of freedom are those 
corresponding to the zero modes of the a fields and the low lying energy levels are 
the eigenvalues of the effective Hamiltonian [21] 

where L 2 is the quadratic Casimir operator with eigenvalues J(J + 1) and 6 c fr is 
the effective moment of inertia (which depends on the volume £). More precisely, 
the lowest energy levels in the isospin J sector are given as®| 

ej(£)=£Ej(£)^g{£)J(J+l), (4.3) 

where g(£) is some effective running coupling. This effective rotator spectrum is 
valid in perturbation theory (PT) up to 2- loop level [221 [23] • The lowest energy 
levels correspond to the totally symmetric tensor states, which in our notation are 
the J = N states (diagonal scattering or £7(1) sector). 

The J = 1 case corresponds to the mass gap. In this case the PT result is known 
up to 3-loop order [24|: 

£i(£) = 27ra(l + a 2 + 1.19 a 3 ) + 0(a 5 ), (4.4) 

where a is a convenient running coupling defined by 

- + lna = ln^ + r'(l)-l. (4.5) 

a I 

Here the constants on the right hand side of this equation are chosen such that the 
1-loop term in (14. 4p vanishes. For general isospin no direct calculation exists beyond 
1-loop order [22], but using the validity of the rotator spectrum up to 2- loop order 
[23] we have 

£j(£) = J(J + l)7ra(l + a 2 ) + 0(a 4 ). (4.6) 
Using the numerical results in Table [lj we established that 

e 2 (£) = 6vra(l + a 2 ) + 0(a 5 ), (4.7) 



8 Note that Ej{£) — Ej(£) — Eq(£), i.e. the energy eigenvalues are measured from the lowest 
energy state (vacuum). 
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£ 




Figure 3: The UV spectrum of the 0(3) model as function of the running coupling a, corre- 
sponding to the range 0.000001 < I < 0.1. The two low energy data sets contain the J = 1 and 
J = 2 rotator energy levels and the corresponding solid lines are the PT predictions (|4.4p and 
(|4.7[) . respectively. Far above the rotator spectrum energies of the N = 2, J = 1 excited states are 
shown. Here the solid line is the constant Air. 



i.e. the 3-loop term in (14.71) vanishes or is very small. 
The excited states above the rotator spectrum start at 

e*(£)^4n, (4.8) 

corresponding to two massless particles of momentum 

The UV spectrum of the a model, obtained by numerically solving the NLIE 
equations down to £ = 10~ 6 agrees very well with the above rotator picture. This is 
shown in Figure [3j 



The Bajnok-Janik formula 

Recently |4] a generalization of Luscher's well-known result [3] about the finite size 
corrections of particle masses has been suggested. According to this conjecture, in 
2-dimensional relativistic integrable modelsj one first defines the auxiliary function 

m) = 7T / dye- 6c ° ah vA(y + iq\e), (4.9) 

9 For the modifications necessary for the nonrelativistic model relevant in the AdS/CFT problem 
see g]. 
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where 9 = (9i, 9 2 , ■ ■ ■ , On) are the physical rapidities of an iV-particle state and 
A(0|£) is the eigenvalue of the transfer matrix corresponding to the state in question. 

Then the leading (for large £) finite size correction to the particle rapidieties are 
given by the solution of the modified Bethe-Yang equations 

e it** 9i A{eM) |l +^| = -1, j = 1,2,..., N. (4.10) 

Finally the energy formula, which includes the leading FS corrections is given by 
E = J2cosh9 j -— dy coshye- ecoshy A(y + iq\9). (4.11) 

j=l ^-0° 

Here the particle rapidities 9j are solutions of the modified BY equations (14.101) . 

Since the above conjecture has been questioned for states with non-diagonal 
scattering in the case of the 0(4) model [2], we decided to numerically investigate 
the validity of it in our model. The results are shown in Tables [3J16J In all cases 
we studied (2-particle states with isospin 0,1,2 and a 3-particle state with isospin 1) 
(I4.10p and (14. lip seem to correctly give the leading FS corrections for the particle 
rapidities and the energy of the state, respectively. 



t 


H (0) 




H 


£(0) 


em 


E 


1 


1.322083 


1.316882 


1.317113 


4.017806 


3.965831 


3.969807 


2 


0.912509 


0.911777 


0.911788 


2.892078 


2.881859 


2.882081 


3 


0.705671 


0.705531 


0.705532 


2.518982 


2.516186 


2.516206 


4 


0.576883 


0.576852 


0.576852 


2.342126 


2.341269 


2.341271 



Table 3: Test of the Bajnok-Janik conjecture for zero momentum N = 2, J = 2 states. is 
the rapidity from the BY equation, ifW from the modified BY equation (|4.10[) and H is the exact 
rapidity. The energy values E^°> , and E are defined analogously. 
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t 


H (0) 




H 


E (0) 




E 


1 


2.424832 


2.426304 


2.426517 


11.388820 


11.579614 


11.606615 


2 


1.756669 


1.757402 


1.757428 


5.965728 


6.006502 


6.008211 


3 


1.392858 


1.393083 


1.393086 


4.274704 


4.285133 


4.285278 


4 


1.155252 


1.155316 


1.155317 


3.489801 


3.492757 


3.492771 



Table 4: Test of the Bajnok-Janik conjecture for zero momentum N = 2, J = 1 states. H^°> is 
the rapidity from the BY equation, ffW from the modified BY equation (|4.10|) and if is the exact 
rapidity. The energy values E(°\ and E are defined analogously. 



i 


H (0) 




H 


£(0) 




E 


2 


1.444416 


1.448838 


1.448743 


4.475261 


4.464046 


4.465255 


3 


1.055912 


1.057006 


1.056993 


3.222472 


3.212762 


3.212965 


4 


0.818360 


0.818622 


0.818621 


2.707934 


2.703854 


2.703885 


5 


0.661971 


0.662035 


0.662035 


2.454443 


2.452936 


2.452940 



Table 5: Test of the Bajnok-Janik conjecture for zero momentum N = 2, J = states. H^°' is 
the rapidity from the BY equation, ffW from the modified BY equation (|4.10p and if is the exact 
rapidity. The energy values E(°\ E^ and E arc defined analogously. 

A Construction of the T-system and TQ-system 

In this appendix we construct T-system elements satisfying ( 12.71) . (12.81) for k > 1, 
(123)1) and (12TT0D for k > 2 and the modified relations (12T2TD and (T2~^2D . assuming 
that a solution of the 0(3) Y-system is found satisfying (11.31) and the Y-functions 
have roots corresponding to the sets (12.161) . 

We start by constructing the functions T 2 and T 3 satisfying 

f k (6 + tq)f k (6-iq) = Y k (6) (A.l) 

for k = 2, k = 3 with the additional assumption that they are bounded for large 9 
and their roots are the elements of the sets t 2 ? ^3- The solution of this problem is 
the fundamental problem in the theory of the TBA integral equations and is given 
explicitly by (k = 2, 3) 

(f) — y-( Q ) \ ( 1 r 00 H?/ 1 
^)^{rJ^7^e^) 1IlYk(u) l' (A ' 2) 
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t 


H (0) 




H 


E (0) 




E 


1 


2.441754 


2.443044 


2.443196 


12.580190 


12.762041 


12.786246 


2 


1.786837 


1.787192 


1.787205 


7.138026 


7.176748 


7.178369 


3 


1.425225 


1.425320 


1.425321 


5.399247 


5.409763 


5.409915 


4 


1.185639 


1.185665 


1.185665 


4.578328 


4.581482 


4.581498 



Table 6: Test of the Bajnok-Janik conjecture for zero momentum N = 3, J = 1 states. H^ ' is 
the rapidity from the BY equation, ffW from the modified BY equation (|4.10|) and H is the exact 
rapidity. The energy values E(°\ and E are defined analogously. 



where t k °^ are the elements of t k . We now define T 4 by. 

fA6) = (A.3) 
T 2 (6) 

Using the Y-system equations and (12.161) it is easy to show that the roots of T4 

are the set £4 and it satisfies (1A.1I) with k = 4. Similarly we construct T 5 ,T e , . . . 
recursively. We also define 

U9) = !^ (A.4) 
and see that its roots are in ti and it satisfies 

Ti(0 + iq) Ti(0 - ig) = ^(0) y o (0). (A.5) 
So far we have constructed Ti, T 2 , . . . satisfying 

f k {e + iq)f k (e-iq) = Y k (9) [l + 6 kl y {6)] (A.6) 



for > 1, 

f k „ l (e)f k+1 (e)=y k (9) (A.7) 
for k > 2 and the T-system equations 

T\(0 + ig) T fe (0 - zg) = 1 + f k+1 (6) f^ff) (A.8) 

for k > 2. The T-system functions can be completed by Tq by requiring (1A.8I) to 
hold for = 1 also. Note that T does not have any roots in the physical strip. 
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The last step is to perform a gauge transformation with g(9) = ip (6 + iq). This 
gives 

<f>(6) = <p {e)<p {e + 2iq) (A.9) 

and 

T k (6) = MO + i(k + l)q) MS - Kk + l)q) f k (0). (A.10) 

Note that T k has the same set of roots in the physical strip (t^) as and it is 
a smooth deformation of the correspnding T-function in the large I Bethe Ansatz 
solution. 

Having constructed the T-system elements we now turn to the TQ-relations 
(12.23p . It is possible to show that they are equivalent to the single second order 
difference equation [25] 



- 2iq) + <j>{6 - 2iq) Q(6 + 2iq) = A{6) Q(0), (A.ll) 

where 



...~)T (e-3iq) + <l>(6-2iq)T 2 (6-iq) 
{ ' T x {9-2iq) ' 1 ' ' 

In the large i limit this can be further simplified: 

A(6) = T 1 (6). (A.13) 

Throughout this paper our main assumption was that the exact solution is a 
smooth deformation of the one obtained from the Bethe Ansatz at large i. Here we 
make this statement more precise. Since the exact TBA equations differ from the 
ones valid in the large i limit only by terms related to Y , which is exponentially 
small in the physical strip, we obviously have 

Y k (6)~Y* A (6), |Im0|<|, A; = 1,2,..., (A.14) 

where ~ here means up to exponentially small corrections (in £) and we introduced 
the superscript BA for objects of the Bethe Ansatz solution. Next, from the set of 
Y-system equations (\1.3\i we see that the deformation is actually exponentially small 
in the larger domain 

n(fl)~lf(fl), M<y> k = l,2,.... (A.15) 
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Inspecting the relation between the T-system and Y-system functions we can see 
that 

T k (9)~T* A (9), |l m ^|<^±l^, k = 1,2, (A.16) 

Finally from the TQ-relations (12.231) we see that it is natural to assume that 

Q(9) ~ Q BA (6) = Qq(9), |Im0|>O. (A.17) 

Thus Q{9) is close to the polynomial Qo(0) in the upper half plane (and its complex 
conjugate Q(9) in the lower half plane). 

In the large I limit it can be shown [2n] that the two linearly independent solu- 
tions of (1A.11I) are Qo(9) (which is a polynomial of degree M) and Rq{9), an other 
polynomial of degree 2N + 1 — M > M. We will assume a similar polynomial 
behaviour of Q(9) in the upper half plane for large 9 thus we can characterize the 
smoothly deformed solution Q(9) uniquely by the requirement that (in the upper 
half plane) 

Q(9) ~ 9 M , (A. 18) 

asymptotically for large 9. 

B NLIE in Fourier space 

In this Appendix we derive the NLIE in the special case where the functions T\ and 
T 2 have only real roots in the physical strip and also the roots of Q(9) satisfy 

\lmu J \< 1 + ^, Q( Uj ) = 0, 3 = l,2,...M. (B.l) 

The set of real Ti roots will be denoted by {rj}. This is the union of the set of 
physical rapidities with the set of type I holes. This important special case covers 
all 2-particle states discussed in this paper (at least for large enough volume) and 
many other multi-particle states of interest. The derivation in the most general case 
is more complicated, but goes essentially along the same lines. 

We start with some definitions. To any bounded meromorphic function ip{z) we 
associate the Fourier transform of its logarithmic derivative along a line parallel to 
the real axis: 

Mz) $(k,a)= [°° dxj kx ^( X + ia } . (B.2) 
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For the complex conjugate we have: 



if)(z) i>(k, a) = i>*{— k, — a). 

The function ip(z) has roots at and poles at P v . We have 



(B-3) 



e- ka ^(k } a) - e- k(3 4>(k } (3) = 2ml ^ e ikR » - ^ v 

^a<Imfl M </3 a<lmP l/ </3 

In the limit (3 — > +oo we have 



ikP v 



(B.4) 



e~ ka ^(k } a) = 2m <{ ^ ~ (> 

Imi?, (J >o ImP„>a 



ifcP„ 



(B.5) 



We start from 



and 



h(9) 
Bi{B) 



which we recall here for k = 1: 

Q(9 + 3iq)T 1 (8-iq) 
Q(9 - 3iq) <f){9 + iq) ' 
Q(6 + iq) T 2 (6) 
Q{9 - 3iq) <f>(9 + iq)' 



;b.6) 
;b.7) 



b 1 {9)b 1 {9) = £,(9) = [l+y 1 (9)}[l + e~ icoshe y 1 (9)}, (B.8) 
B 1 (9 + iq)B 1 (9-iq) = £ 2 (9) = l + y 2 (9). (B.9) 



We now define 



n(Jfe) = 2vrz ^ e ifc ^ , r 2 ( jfe) = 2vri ^ e ifc ^ 



(B.10) 



Some further definitions: 



r(k)=e-" k f 1 (k,q)-f 2 (k,q), 



(B.11) 



= 0), £ 2 (A;) = £ 2 (fc,0), &(*) = 6i(fc, 7 ), B(k) = B,(k, 7 ). (B.12) 
^From the ratio of (1B.6P and (IB. 71) we have 



6(fc) = B{k) + f^fc, T - q) + Q(k, 7 + 3g) - T 2 (fc, 7 ) - Q(fc, 7 + g). (B.13) 
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Since Q(9) is analytic in the upper half plane and has no roots in [7 + q, 00] we 
conclude using (1B.5P that 

Q(k, 1 + q) = Q(k, 1 + 3q) = } k>0. (B.14) 

Furthermore using (1B.4j) we see that 



f l (k, 1 -q)=e k ^f 1 (k,q)+e^ k T 1 (k), f 2 (k, 1 )=S^f 2 (k,q) (B.15) 

and thus we can write 

b{k) = B(k) + e^~ q)k [r{k) + T^k)] , k>0. (B.16) 

We now use the facts that 
T\{6) has no roots in [— q, — q + 7], 
Q{9) has no roots in [3q, 3q + 7], 
(p(9) has no roots in [q, q + 7], 

has no roots in [— 3q, —3q + 7], 
and conclude from (IB. 61) using (IB. 41) that 

6i(fc, 0) = e - fc7 6(A:) (B.17) 

and using (1B.8|) that 

&i(/c, 0) + bl(-k, 0) = li(k) = e- k ~ f b{k) + e k W{-k). (B.18) 

This can be used to write a relation also for negative k: 

b(k) = e fc7 £i (k) - e 2k ^B*{-k) - e k ^ +q) [r*{-k) - n (k)] , k < 0. (B. 19) 

Similarly using that 
T 2 {9) has no roots in [7, q], 
Q{9) has no roots in [7 + q, 2q], 
(f)(9) has no roots in [7 + q, 2q], 
Q{9) has no roots in [7 — 3q, —2q], 
we conclude from flB.71) using flB.41) that 

B 1 (k,q) = e k{q - l) B(k) (B.20) 
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and also using (1B.9|) we have 

£ 2 (k) = e k ^B(k) + e k ^~ g) B*(-k). (B.21) 
We now go back to the original relations (12.101) and (I2.22p and write 

h{k) = f%(k, q) + f^k, -q) - 4>(k, q) - }(k, -q) (B.22) 

and 

£ 2 (k) = f 2 (k, q) + f 2 {k, -q) - 4>{k, 2q) - }{k, -2q). (B.23) 
Since (f)(9) has no roots in [— 2q, —q] and (f)(9) has no roots in [q, oo] we see that 

e 2qk }(k, -2q) = e qk }(k, -q) (B.24) 

and 

4>(k, 2q) = 4>(k, q) = 0, k > 0. (B.25) 

Further we note that 

e kq f % (k, -q) - e~ qk f % (k, q) = r t (k), i = 1, 2. (B.26) 
Combining (IB. 22[) and (IB.23[) we thus have 



e- qk h(k) - £ 2 (k) = r(Jfe) (1 + e~ 2qk ) + e~ 2qk n(k) - e~ qk r 2 (k), k > 0, (B.27) 
and, after some algebra and using (lB.2ip 

r(k) = g(k){l 1 (k) - e k{2q ~^B(k) - e k ~<B*(-k) 

(B.28) 

+ T2 (A;)-e-* fc ri(A;)}, k > 0, 

where 

**> = 2^ < R29) 

^From this we get for negative k, after complex conjugation: 
r*(-k) = n(Jfc) + g(k){! 1 (k) - e~ k W(k) - e k{ ' ( - 2q) B*(-k) 



(B.30) 

-T 2 (k)-e- qk Tx(k)j, k<0. 



28 



Finally we use (IB. 28[) in (IB. 16j) and similarly (IB. 301) in (IB. 19j) and get for positive 
k: 



b(k) = g(k)\e~ qk B(k) - S 2 ^B*(-k) 

1 , (B.31) 

+ e^-^l^k) + e k ^ q) T 2 (k) + e k ^n(k)\, k > 

and for negative k: 

~b(k) = g{k)\e qk B{k) -e k ^ +q ^B*{-k) 

, (B.32) 

+ e^-^e^k) + e k ^ +q) r 2 (k) + e^nik) \, k < 0. 



Combining the last two equations we now write the NLIE in Fourier space as 

b(k) = g(k)\e- qW B(k) - e 2k ^- q ^B*(-k) 

1 , (B.33) 

+ e fc(7 - g) £i(fc) + e fc7 - <?|fc|) r 2 (fc) + e fc7 n(fc) \. 



This is the Fourier space version of the first equation of the NLIE (12.301) in our 
special case. 
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